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ABSTRACT 

We present a systematic analysis of the X-ray spectral properties of a sample of 
22 "narrow-line" Seyfert 1 galaxies for which data are available from the ASCA public 
archive. Many of these sources, which were selected on the basis of their relatively 
narrow H/3 line width (FWHM <2000 km/s), show significant spectral complexity 
in the X-ray band. Their measured hard power-law continua have photon indices 
spanning the range 1.6 — 2.5 with a mean of 2.1, which is only slightly steeper than 
the norm for "broad-line" Seyfert Is. All but four of the sources exhibit a soft excess, 
which can be modelled as blackbody emission (Tbb w 100 — 300 eV) superposed on 
the underlying power-law. This soft component is often so strong that, even in the 
relatively hard bandpass of ASCA, it contains a significant fraction, if not the bulk, of 
the X-ray luminosity, apparently ruling out models in which the soft excess is produced 
entirely through reprocessing of the hard continuum. 

Most notably, 6 of the 22 objects show evidence for a broad absorption feature 
centred in the energy range 1.1-1.4 keV , which could be the signature of resonance 
absorption in highly ionized material. A further 3 sources exhibit "warm absorption" 
edges in the 0.7-0.9 keV bandpass. Remarkably, all 9 "absorbed" sources have H/3 line 
widths below 1000 km/s, which is less than the median value for the sample taken as 
a whole. This tendency for very narrow line widths to correlate with the presence of 
ionized absorption features in the soft X-ray spectra of NLSls, if confirmed in larger 
samples, may provide a further clue in the puzzle of active galactic nuclei. 
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1 INTRODUCTION 

Significant hard X-ray luminosity is now regarded as one of 
the defining characteristics of Seyfert 1 galaxies. X-ray spec- 
tral and variability observations establish that the emission 
is produced in the innermost regions of the Seyfert nucleus 
and provide a means of investigating both the nuclear ge- 
ometry and the processes accompanying accretion onto a 
supermassive black-hole (e.g. Mushotzky, Done & Pounds 
1993). 

The hard X-ray spectra of normal Seyfert Is can, to first 
order, be described in terms of a power-law continuum with a 
photon index of typically F = 1.9 (Nandra & Pounds 1994). 
In many sources additional spectral features are observed, 
such as a spectral turn-up above 10 keV and strong iron Ka 
fiuorescence in the 6-7 keV band, which can be attributed 
to Compton reflection of the continuum by optically thick 
matter, possibly in the form of an accretion disk (e.g. George 
& Fabian 1991; Matt, Perola & Piro 1991). The presence of 



a soft X-ray excess may in turn represent the high energy 
tail of the thermal emission generated in the inner regions of 
the disc (e.g. Czerny & Elvis 1987). Recent interest has also 
focussed on the detection of "warm absorption" features in 
the soft X-ray spectra of Seyfert 1 galaxies which are an 
imprint of highly photoionized material lying in the line of 
sight to the nuclear source (Reynolds 1997; George et al. 
1998). 

Narrow- line Seyfert Is (NLSls) represent a subclass of 
objects occupying one extreme of the measured range of op- 
tical line widths, i.e. H/3 FWHM < 2000 km/s (Osterbrock 
& Pogge 1985). Observations with ROSAThme shown that 
NLSls frequently have stronger soft excesses and increased 
variability in the 0.1-2 keV band than more normal, "broad- 
line" Seyfert Is (hereafter BLSls; Boiler, Brandt & Fink 
1996). A possible explanation of the extreme properties of 
NLSls is that these are active galactic nuclei (AGN) con- 
taining black holes of relatively modest mass which, never- 
theless, are accreting at a high rate (Pounds, Done & Os- 
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borne 1995). In this model, the lower mass black hole gives 
rise to smaller Keplerian velocities of the broad-line region 
clouds while the high accretion rate yields comparable lu- 
minosity to that of larger mass black holes and leads to 
increased disc emission and an enhanced soft excess (Ross, 
Fabian & Mineshige 1992). In some NLSls there is tenta- 
tive evidence for reflection from a highly ionized disc (e.g. 
Comastri et al. 1998), which fits naturally with the high 
accretion rate model, as the disc surface is expected to be 
highly ionized (Matt, Fabian & Ross 1993). Recent studies 
suggest that NLSls, in addition to their extreme soft X-ray 
spectra, also exhibit steeper intrinsic hard X-ray continua 
than their broad-line counterparts (e.g. Brandt, Mathur & 
Elvis 1997), an effect which could arise, for example, if the 
strong soft excess effectively cools the accretion disk corona 
in which the underlying power-law is formed (Pounds, Done 
& Osborne 1995; Maraschi & Haardt 1997) 

Our objective in this paper is to investigate the X-ray 
spectral properties of a sample of NLSls. Specifically we are 
interested in whether there are any further X-ray spectral 
characteristics of NLSls, over and above those noted ear- 
lier, which might be described as distinctive of this class 
of object. In practical terms we have carried out a system- 
atic analysis of 24 observations of 22 NLSls taken from the 
ASCA public data archive. The remainder of this paper is 
organised as follows. In the next section we define the source 
sample and detail the observations and data reduction tech- 
niques. In Section 3 we go on to describe the results of our 
X-ray spectral analysis and in Section 4 present a discus- 
sion of these results. Finally we summarise our conclusions 
in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

We have selected 22 objects, identified in the literature as 
Seyfert 1 galaxies or quasars with H/3 FWHM <2000 km/s, 
for which ASCA data were available at the time of the analy- 
sis (late 1998). The sample, given in Table 1, is not complete 
in any sense except that it consists of all suitable spectra of 
NLSls in the ASCA public archive. The redshift distribu- 
tion of the sample is shown in Fig. 1. The NLSl galaxy 
Mkn 957 was also observed by ASCA but these data have 
not been analysed due to the very short exposure time of 
the observation (~3 ks). In addition data on Kaz 163 have 
been excluded; Kaz 163 is in the same GIS field as Mkn 507 
but its image falls at the extreme edge of the SIS detector, 
making the data unsuitable for detailed spectral analysis. 

ASCA carries four instruments which are operated si- 
multaneously namely two solid-state imaging spectrometers 
(SIS-0 and SIS-1; Burke et al. 1991) and two gas imaging 
spectrometers (GIS-2 and GIS-3; Kohmura et al. 1993). In 
the current analysis, standard data screening criteria were 
applied to all the data. Table 1 lists the resulting expo- 
sure times (specifically for SIS-0, although data from all 
four ASCA instruments are utilised here). Since we are con- 
cerned only with the spectral properties of the sources, time- 
averaged spectra have been extracted from each observation. 
Counts were accumulated, where possible, in circular aper- 
tures of 3-4 arcmin radius centred on the source. The back- 
ground was estimated using source free regions from the 
same observation and, in the case of the GIS observations. 
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Figure 1. The redshift distribution of the NLSl sample. 



at a similar off-axis angle to the source. After background 
subtraction the pulse height spectra were binned to give at 
least 20 counts per spectral channel. In the subsequent anal- 
ysis the SIS data were restricted to the 0.6-10 keV energy 
range and GIS data to the 0.8-10 keV band so as to reduce 
the possible impact of calibration uncertainties. 

As a check of the possible impact of spectral variability 
on the time-averaged source spectra, we extracted the 2-10 
keV and 0.5-2.0 keV light curves of each source (SIS-0 and 
SIS-1 co-added) and binned these on the satellite orbital 
timescale (5760 s) . We then used a test to check the con- 
stancy of the spectral hardness (defined as ratio of the 2-10 
keV count rate to the 0.5-2.0 keV count rate) during the ob- 
servation. Five out of 22 sources gave unacceptable values of 
at the 95% level (NAB 10205-^024, RX J0439-45, NGC 
4051, Mkn 478, Ark 564) indicative of underlying spectral 
variability. However, with the exception of NGC 4051 there 
was little evidence that these variations corresponded to any 
clear spectral trends or to a systematic effect (e.g. a spectral 
hardness versus flux correlation). We conclude that the im- 
plicit assumption of the analysis which follows, namely that 
the observation-averaged spectra provide a good represen- 
tation of the source properties, is a reasonable one. 



3 X-RAY SPECTRAL ANALYSIS 

Spectral fltting of the binned source spectra was carried 
out using the xspec vlO.O software package. In general, 
data from all four detectors have been fitted simultane- 
ously, but with the relative normalisations of the spectra 
free to vary between the four instruments. The quoted er- 
rors on the derived best-fitting model parameters correspond 
to a 90% confidence level for one interesting parameter (i.e. 
a Ax^ = 2.7 criterion) unless otherwise stated. Values of 
Hq = 50 km/s/Mpc and go = 0.5 are assumed throughout 
this section and fit parameters (specifically the emission and 
absorption line energies) are quoted for the rest frame of the 
source. 
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Table 1. The ASCA NLSl sample. The columns contain the following information: (1) The source name; (2) The sequence number 
of the observation; (3) & (4) The observation pointing position; (5) The exposure time for SIS-0 (ks); (6) The line-of-sight Galactic 
hydrogen column density from Dickey & Lockman (1990) (in units of lO'^" cm~^); (7) The source redshift; (8) The H/3 line width (km/s); 
(9) References to the information in columns (7) & (8). 
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REFERENCES; (1) Wang ct al. (1996). (2) Brandt (1996). (3) Grupe (1996). (4) Zheng & O'Brien (1990). 



3.1 Initial spectral fits 

Broadly, our strategy has been to apply relatively simple 
spectral models, so as to characterize the properties of the 
whole sample of sources. Thus initially a model comprising 
just a power-law continuum with absorption due to the line- 
of-sight Galactic Nh (see Table 1) was employed as a way of 
illustrating the variety of X-ray spectra exhibited by NLSl 
galaxies. Fig. 2 shows the residuals from these fits for all 
24 observations, which in most instances have unacceptable 
minimum values. Cursory inspection of Fig. 2 shows that 
soft excesses are often present below ~2 keV and that broad 
absorption and/or emission complexes are also evident in a 
number of the sources. In general, an iron Ka emission-line 
is difficult to pick out (by eye), although NGC 4051 provides 
the exception. In one source, Mkn 507 the largest residual 
occurs in the lowest spectral channel suggesting the possi- 
bility of cold absorption over and above the Galactic column 
density. In this case we obtain an improved fit (A^^ ~ 12) 
when the Nh is allowed to be a free parameter, the specific 
value derived being Nh = 3.2 ± 1.5 x lO'^^ cm"'^ (see Iwa- 
sawa, Brandt & Fabian 1998). In the rest of this paper we 
assume that Mkn 507 is the only NLSl in the sample with 
excess cold absorption. 

3.2 The 2-10 keV spectral form 

In order to characterize better the hard X-ray spectrum of 
each NLSl, we next excluded all the spectral data below 



2 keV, where evidently soft excesses and warm absorbers 
add to the spectral complexity. The 2-10 keV spectra were 
again fitted with a simple power law plus Galactic absorp- 
tion model but with the addition of a narrow Gaussian emis- 
sion line [a = 0.01 keV) to represent iron Ka emission in 
the 6.4-7.0 keV bandpass. Table 2 summarizes the results. 
Acceptable fits were found for most sources, with only three 
instances of (reduced) > l-l- We note that the derived 
photon indices have a mean value of 2.1 with a standard 
deviation of 0.3. 

For nine objects the inclusion of the iron-line compo- 
nent is merited in terms of the resulting improvement in 
{i.e. Ax^ > 6.2 implying better than 95% confidence in 
the F-test for two additional parameters). Also two further 
sources show at least a marginal improvement in (see 
Table 2). There is some evidence for lines originating from 
fairly highly ionized iron species (i.e. >Fexx and above) in 
I Zw 1, Ton S180, PKS 0558-504, PG 1244+026 and Ark 564. 
The measured line equivalent widths range from 100-600 eV, 
albeit with large uncertainties. Similarly the constraints on 
the equivalent width of a (neutral) iron line in those sources 
lacking a significant line detection are generally rather weak. 
In general the signal/noise ratio was too poor to meaning- 
fully constrain the intrinsic line widths. However, the second 
observation of NGC 4051 provides the exception since in this 
case there is evidence for a broadened iron Ka feature. The 
best fit parameters for the line are E= 6.29 ± 0.09 keV, in- 
trinsic width a = Q.mto'S, keV and EW= 2AQtfo eV. The 
improvement in the fit upon adding this broad feature is 
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A.x^ — 76.8. An attempt has been made to fit the hne with 
the relativistic "diskline" profile of Fabian et al. (1989); we 
find a mildly ionized line (E= 6.6lg2 keV), with a near 
to face-on inclination (i < 17°) gives the best fit to the 
data (assuming values for Ri„, Rout and emissivity index 
of 10GM/c^ lOOGM/c^ and -2, respectively). Note that 
Guainazzi et al. (1996) discuss the the spectral fitting of 
this ASCA observation of NGC 4051 in some detail. 



3.3 The soft X-ray spectra 

Extrapolation of the best-fit 2-10 keV spectrum (as defined 
in Table 2) down to 0.6 keV in most cases results in a poor fit 
of the soft X-ray spectrum, with the most common residual 
feature being an excess of soft flux. In the spectral fitting 
we have attempted to match this soft excess with an ad- 
ditional continuum component. Specifically we use a single 
blackbody component, although a second power law often 
provides an equally good fit. (In the latter case the second 
power-law is typically steeper than the first by AF ~ 0.5 
with a break energy in the range 1-2 keV). Table 3 sum- 
marises the results of fitting a power law plus blackbody 
model (note that from here on the iron line parameters are 
frozen at the values obtained in the earlier 2-10 keV fits). 
In all but four objects (I Zw 1, PG 1543-^489, Mkn 507 
and IRAS 20181-224) a soft excess component provides a 
significant improvement in the fit, demonstrating that soft 
excesses are a very common feature in NLSls. The underly- 
ing power-law photon indices given in Table 3 have a mean 
of 2.12 and a standard deviation of 0.26. The fact the mean 
is very similar to that obtained earlier for 2-10 keV fits 
demonstrates that the latter fits are relatively immune to 
the presence of the soft excess. 

Nine objects show signs of additional spectral complex- 
ity below 2 keV even after the inclusion of the blackbody 
component in the fit (Table 3). Unfortunately modelling of 
the ASCA spectra in terms of additional soft X-ray features 
is not particularly straight-forward. Radiation damage to 
the CCDs and other factors have meant that there is in- 
creasing uncertainty in the calibration of the detectors be- 
low 1 keV and, in particular, it has been noted that the 
two SIS instruments often give divergent spectra at the low- 
est energies (although these calibration uncertainties only 
dominate over uncertainties in the background subtraction 
for relatively bright sources). Also since most NLSl galaxies 
appear to exhibit a soft X-ray excess it is difficult to dis- 
tinguish subtleties in the form of the soft continuum from 
the eflects of putative broad emission and/or absorption fea- 
tures. 

Recent studies of the X-ray spectra of NLSls (e.g. 
Leighly et al. 1997b; Fiore et al. 1998) have established that 
in addition to classical "warm-absorption" features, NLSls 
often show anomalous absorption features in the 1-2 keV 
band. Both types of absorption are spectrally complex and 
will merit more detailed analysis using the predictions of ap- 
propriate photoionization codes, once high sensitivity X-ray 
data with good spectral resolution become available from 
missions such as AXAF, XMM and ASTRO-E (e.g. Nicas- 
tro et al. 1999). However, for our present purpose we have 
taken a very simplistic approach and have attempted to im- 
prove the in the spectral fits for nine sources noted above 



Table 4. X-ray spectral fitting of a Gaussian absorption line. The 
columns give the following information: (1) The source name; (2) 
The best fit line energy (keV); (3) The line equivalent width (eV); 
(4) The intrinsic line width (keV); (5) The improvement in the fit 
when the absorption feature is included (compared to Table 3). 



Name 
(1) 

Ton S180 
IH 0707-495 
NGC 4051 
PG 1244+026 
IRAS 13224-38 
PG 1404+226 
IRAS 17020+45 
IRAS 20181-22 
Ark 564 



Energy EW 
(2) (3) 



1.22"* 



1.31"* 



,+0.07 
-0.04 

1.08 ± 0.03 
0.83 ± 0.03 

1+0.07 
-0.03 

1.10 + 0.07 
1.16 ± 0.05 

0.90 + 0.03 
1.38 + 0.03 



-95 + 16 
-112^ 
-33^ 



,+25 
-38 
3+4 
-21 

-108 + 50 

H20 
-160 

-100 + 30 

^^^-35 

_42+6l 
^^-24 
-+24 
-10 



(4) 

0.31 ± 0.04 
0.14 + 0.07 
0.15 



-275"* 



-55^ 



+0.15 
-0.02 

0.24 + 0.08 

^•^^-0.03 

0.09 + 0.02 
0.12 + 0.05 
< 1.6 



0.19 



+0.11 



(5) 

20 
32 
43 
24 
87 
31 
64 
8 
132 



by including just a single absorption feature (in the form of 
a broad Gaussian absorption line) in the spectral model. 

Table 4 lists for each of the nine sources the line energy, 
the equivalent width and the intrinsic line width obtained 
when such an absorption feature is added to the power-law 
plus blackbody continuum model. The line energies are not 
consistent with a single value but a bifurcation is suggested, 
namely absorption either in the 0.7-0.9 keV range or in the 
1.1-1.4 keV range. Absorption features in the former range 
are usually interpreted as due to Ovil and Ovill edges in- 
dicative of the presence of ionized material along the line 
of sight. The possible origin of the "anomalous" absorption 
features observed at ~ 1.2 keV have recently been discussed 
by Leighly et al. (1997b) and Fiore et al. (1998). 

As noted above the fitting of a Gaussian absorption fea- 
ture is necessarily an over simplification of the true picture. 
In order to investigate the individual sources in somewhat 
more detail we have refitted the spectra of the nine objects 
listed in Table 4 including a variety of additional absorp- 
tion and emission components. For example, we tested for 
either one or two absorption edges and also for either a single 
Gaussian emission line or a MEKAL-type optically thin ther- 
mal spectrum (Kaastra & Mewe 1993). The outcome was 
that in five of the objects an absorption component gave 
a significantly better fit than an emission component. (Of 
course the modelling of the underlying soft excess changes 
substantially between these two cases so as to maintain the 
match to the observed spectrum in the 0.6-2 keV bandpass). 
Unfortunately the situation is more ambiguous in the re- 
maining cases since absorption and emission models can be 
constructed which give fairly comparable fits to the ASCA 
data. Further comments on the spectral fitting of the indi- 
vidual nine sources are given in the Appendix. 



4 DISCUSSION 

4.1 The X-ray continua of NLSls 

Using spectral measurements compiled from the published 
literature Brandt et al. (1997) identify a correlation between 
the hard power-law spectral index and the width of the opti- 
cal H/3 line in a sample of Seyfert 1 galaxies, in the sense that 
NLSls tend to have somewhat steeper hard X-ray spectra 
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Table 2. The results of the 2—10 keV spectral fitting. The columns give the following information: (1) The source name; (2) The 
powor-law photon index; (3) The 2-10 keV flux (10"^'^ erg cm~^ s~l); (4) The 2-10 keV luminosity (10*^ erg s"^); (5) The iron line 
energy (keV); (6) The iron-line equivalent width (eV); (7) The improvement in the fit when the iron line is included; (8) The reduced 
and number of degrees of freedom for the best fitting model. 

Ta-io Flux2_io L2-10 Ef-e EWfe Ax|,^ X^/dof 

(1) (2) (3) (4) (5) (6) (7) (8) 
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5.9 


Ton SI 80 


2.39 


± 





.05 


45.3 


8.3 


PHL 1092 


1.67 


± 


0, 


.30 


3.6 


32.0 
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IRAS 17020+454 


2.20 


± 


0, 


,06 


78.4 


12.6 


Mkn 507 


1.61 


± 


0, 


,30 


6.7 


9.3 


IRAS 20181-224 


2.33 


± 


0, 


,12 


10.4 


17.5 


Ark 564 


2.50 


± 


0, 


,03 


205.7 


5.2 



6.41 ± 0.08 


193 + 78 


16.3 


1.00/458 


6.73+O-^f, 

— 0.30 


450+^5? 

— 225 


11.0 


1.04/262 


6.65 ± 0.10 


180 + 90 
< 490 


10.9 


1.05/591 
1.18/132 


6.5 ±0.2 


670 + 300 


8 


0.86/87 


_ 


< 122 


_ 


0.85/327 




< 220 




1.09/149 


6.65+0-35 

— O.zO 


65 + 40 


6.5 


0.98/820 




< 375 




1.04/114 




< 490 




1.08/98 




< 1560 




0.97/32 




< 380 




0.93/122 


6.43 ± 0.04 


145 + 40 


35.6 


1.00/905 


6.37 + 0.03 


127 + 25 


66.5 


1.06/1449 


f. 00+0.12 


260 + 130 


10.7 


0.94/342 


7.0 + 0.1 


460 + 210 


11 


0.99/255 




< 290 




1.12/182 




< 1220 




0.84/57 




200 + 135 


5.7 


1.02/216 




< 350 




0.87/76 


6.5 + 0.5 


95 + 70 
< 715 


5 


0.87/577 
1.02/50 




< 120 




1.14/185 


7.0 + 0.3 


123 + 50 


15 


0.99/1171 



Table 3. The full 0.6—10 keV spectral fits. The columns give the following information: (1) The source name; (2) The power-law photon 
index; (3) The blackbody temperature (eV); (4) The 0.6-10 keV flux (10"^'^ erg cm~^ s"^); (5) The derived unabsorbod luminosity 
(10*3 erg s~l); (6) The ratio of the 0.6-10 keV luminosities of the blackbody and power-law components; (7) The improvement in the 
fit when the blackbody component is added; (8) The reduced \^ and number of degrees of freedom for the best fit model. 

Name V kT FluxQ.e-io L0.6-10 Lftb/Lpi x5/dof 

(1) (2) (3) (4) (5) (6) (7) (8) 



Mkn 335 


1.96 + 0.04 


160 + 14 


191 


5.7 


0.12 


125 


1.01/757 


I Zw 1 


2.30 + 0.03 


150^ 


72 


14.0 


< 0.01 




0.94/594 


Ton S180 


2.36 + 0.05 


190 + 5 


137 


24.7 


0.29 


226 


1.10/988 


PHL 1092 


1.67 + 0.11 


162 + 15 


7.0 


73 


0.80 


112 


1.12/268 


RX J0148-27 


2.15 + 0.15 


1 29+3^ 


18.3 


12.7 


0.15 


12 


0.91/230 


NAB 0205+024 


2.10 + 0.08 


218 + 19 


53 


63 


0.16 


36 


0.98/697 


RX J0439-45 


2.28 + 0.10 


137 + 11 


27.4 


81 


0.58 


152 


1.02/379 


PKS 0558-504 


2.26 + 0.03 


230 + 15 


257 


247 


0.08 


34 


1.03/1217 


IH 0707-495 


2.33 + 0.10 


107 + 4 


39.6 


3.7 


1.21 


429 


1.29/338 


RE J1034+39 (1) 


2.35 + 0.13 


131 + 9 


32.8 


2.7 


0.64 


104 


0.95/298 


RE J1034+39 (2) 


2.45 + 0.20 


142 + 20 


33.6 


2.8 


0.49 


23 


0.78/125 


Mkn 42 




227 + 40 


20.6 


0.53 


0.15 


8 


0.99/271 


NGC 4051 (1) 


1.97 + 0.02 


86 + 5 


561 


0.10 


0.22 


824 


1.04/1205 


NGC 4051 (2) 


1.94 + 0.01 


106 + 5 


402 


0.07 


0.11 


1160 


1.13/1846 


PG 1211+143 


2.07 + 0.06 


110 + 12 


72 


23.0 


0.17 


57 


0.95/713 


PG 1244+026 


2.31 + 0.08 


210 + 10 


75 


8.1 


0.26 


64 


1.17/535 


IRAS 13224-3809 


1.60 + 0.13 


119 + 4 


18 


4.7 


2.1 


537 


1.28/426 


PG 1404+226 


1.87 + 0.20 


115 + 7 


12.9 


9.1 


1.9 


177 


1.20/167 


Mkn 478 


1.96 + 0.05 


89 + 12 


50 


14.5 


0.23 


101 


1.03/495 


PG 1543+489 


2.48 + 0.20 


436 + 436 


10.7 


105 


0.04 


1 


0.91/206 


IRAS 17020+454 


2.21 + 0.05 


260 + 36 


148 


29 


0.05 


9 


1.06/972 


Mkn 507 


1.74 + 0.24 


150-'' 


8.6 


1.5 


< 0.08 




0.96/97 


IRAS 20181-224 


2.33 + 0.11 


258l?2 


19.3 


47.4 


0.08 


3 


1.15/389 


Ark 564 


2.44 + 0.03 


195 + 4 


568 


16.7 


0.21 


558 


1.27/1567 
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Figure 2. The residuals from spectral fits to the 0.6—10.0 keV ASCA spectra using a power-law plus Galactic absorption model. The 
residuals are in terms of sigmas, with error bars of size one. For clarity only the SIS data are shown, and the SIS-0 data are marked by 
squares. 
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Figure 3. The photon index of the hard continuum plotted 
against the H/3 line width for the NLSl sources in the ASCA 
sample. The dotted line represents the mean spectral index of 
the BLSl sample of Nandra et al. (1997). 



PG 1404+226), the blackbody component contains more lu- 
minosity than the power law in the 0.6-10 keV range, i.e., 
Lbb/Lpi > 1 (see Table 3). For these three sources, extrapo- 
lation of the continuum spectrum down to 0.15 keV and up 
to 100 keV still leaves this ratio above unity. The implica- 
tion is that, at least in these NLSls, the soft excess cannot 
solely be due to reprocessing of the hard X-ray continuum. 
(Note for these three objects the soft excess is still dominant 
even after the inclusion of additional spectral features, see 
the Appendix.) 

Two objects with strong soft excesses, namely 
PHL 1092 and IRAS 13224-3809, appear to have partic- 
ularly hard underlying X-ray continua (F = 1.6 — 1.7, see 
Table 3). This could be used as an argument against mod- 
els where the soft photons Compton scatter against the hot 
electrons in an accretion disk corona to produce the under- 
lying X-ray hard spectrum. However, an alternative view 
might be that the relatively flat 2-10 keV spectra are the 
result of strong reflection in these objects. Sensitive broad- 
band X-ray spectral measurements (extending well above 
the Compton peak) would help distinguish between these 
two descriptions. 



than normal broad-line objects. Fig. 3 shows a plot of power- 
law photon index (from Table 3) against the H/3 FWHM line 
width for our source sample. Clearly within the restricted 
range of hne width which defines the NLSl subclass {i.e. H/3 
FWHM < 2000 km s'^ ) there is no hint of a correlation. As 
noted earlier the mean photon index of our sample of NLSls 
is F = 2.12, with a standard deviation of a— 0.26. Nandra 
et al. (1997) tabulate the 3-10 keV spectral slopes for 15 
BLSls based on ASCA measurements. These give an aver- 
age slope of r = 1.9 with a standard deviation of a = 0.17. 
Thus a modest trend in F over a broad range of H/3 line 
width is suggested, consistent with the Brandt et al. (1997) 
result. However, we caution that the spectral steepening is 
a fairly subtle effect and the above comparison may not be 
entirely free of selection bias. For example, one could argue 
that the comparison is of a hard X-ray selected sample with 
one largely selected at soft X-rays (but see the discussion 
in Brandt et al. 1997). In counterbalance is the fact that 
Nandra et al. (1997) include a reflection component in their 
analysis, which can have the effect of increasing the derived 
power-law spectral indices by up to AF ~ 0.1. Thus it is pos- 
sible that magnitude of the spectral steepening trend noted 
above is understated. 

In the present analysis we have fitted the soft excess in 
terms of a single temperature blackbody component, which 
should be a reasonable assumption, given the restricted 
ASCA bandpass, if the emission resembles that expected 
from an accretion disc (e.g. Ross et al. 1992). Indeed, the 
addition of a blackbody component significantly improves 
the fit (at the > 99% confidence level in the F-test) in 18 
out of 22 objects. This is a much higher rate of incidence of 
soft excesses than appears to be the case for BLSl samples 
(e.g. Turner & Pounds 1989; Reynolds 1997). The soft ex- 
cess generally dominates over the hard power law at energies 
^1.5 keV, and presumably peaks in the extreme ultravio- 
let. Even in the limited ASCA bandpass the luminosity of 
the soft excess is comparable to that of the hard power law 
and in three objects (IH 0707-495, IRAS 13224-3809 and 



4.2 Iron Ka emission 

Emission lines at 6-7 keV are detected at a modestly signif- 
icance level in 9 objects with a marginal indication of such 
lines in 2 others. However, as noted earlier the line proper- 
ties (energy, equivalent width, intrinsic width) are in general 
only poorly constrained by the ASCA spectra. A point of in- 
terest, for follow-up when more sensitive X-ray spectra are 
available, is the incidence of emission lines from highly ion- 
ized gas (corresponding to Fexx and above) in NLSls, for 
which the current data provide only a fleeting glimpse. 

4.3 Soft X-ray spectral features 

The residuals to the best fltting continuum models show un- 
usual features below 2 keV in nine objects. In three sources 
(NGC 4051, IRAS 17020-^454 and IRAS 20181-224) we 
are probably detecting Ovil and Ovill edges at 0.74 and 
0.87 keV respectively, indicative of the presence of "warm" 
absorbing gas in our line of sight to the nuclear contin- 
uum source. Of the other six objects, three (IH 0707-495, 
IRAS 13224-3809 and PG 1404-1-226) are most likely af- 
fected by absorption, whereas for the others (Ton S180, 
PG 12444-026 and Ark 564) the interpretation is more am- 
biguous (see the Appendix). However, for the purpose of 
the present discussion we pursue the assumption adopted 
earlier, namely that all six sources exhibit anomalous ab- 
sorption in the range 1.1-1.4 keV with the equivalent width 
of the absorption feature being typically ~ 100 eV and its 
intrinsic width (when modelled as a pure Gaussian) ranging 
from 0.1-0.3 keV (see Table 4 and Fig. 4). 

Leighly et al. (1997b) have recently reported the de- 
tection of the same anomalous absorption features in three 
sources and discussed possible physical origins of the effect. 
Specifically these authors note that in order to identify the 
absorption with either an Ovil and /or an Ovill edge, as 
predicted in the warm absorber scenario, then very signifi- 
cant (i.e. relativistic) outfiow velocities would be required in 
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Energy (keV) 

Figure 4. The ratio of the ASCA data to the best-fit contin- 
uum model for two sources, IRAS 13224—3809 (filled circles) and 
IRAS 17020-1-454 (open squares). The former source exhibits an 
anomalous absorption feature near 1.2 keV and the latter proba- 
bly a warm (OVII) edge near 0.74 keV. Both spectra have been 
shifted into the source rest frame and for clarity only SIS-0 data 
are shown. 



order to achieve the necessary blueshifting of the edge ener- 
gies. An alternative model has been suggested by Nicastro, 
Fiore & Matt (1999) who invoke resonant absorption lines 
in a warm absorber to explain the absorption evident in the 
spectrum of IRAS 13224-3809, in which case there is no re- 
quirement for relativistic outflow. The basis of their model 
is that a steep soft X-ray spectrum can produce a differ- 
ent ionization structure in the warm absorber, with carbon 
and oxygen fully stripped. Absorption is then produced by 
a complex of absorption lines around 1-2 keV, made up of 
mainly iron L resonance lines (see figure 4 of Nicastro et 
al. 1999). The total equivalent widths of these lines can be 
~100 eV, consistent with our results, although of course our 
assumption of a Gaussian absorption profile is in this case 
an extremely crude approximation to the true underlying 
spectral form. 

On the basis of the above arguments, it is at least plau- 
sible that both the 0.7-0.9 keV and the 1.1-1.4 keV ab- 
sorption arises due to the presence of highly ionized gas 
along the line of sight to the nuclear X-ray source. Interest- 
ingly the two different types of absorbers (described earlier 
as "warm" and "anomalous") occupy very different regions 
of parameter space in Fig. 5, which shows the ratio of the 
luminosities of the soft excess to the hard power law compo- 
nents (Ltb/Lpi) versus the H/3 line width. The sources with 
anomalous absorption clearly tend to have relatively high 
Lbb/Lpi values whereas the three sources exhibiting "nor- 
mal" warm absorption are found at the low end this range. 
This is circumstantial evidence in support of the Nicastro 
et al. (1999) model, in that they require a strong contribu- 
tion from the soft excess in order to produce the necessary 
conditions for the dominance of resonance over pure photo- 
electric absorption. This point is given further weight if we 
also consider BLSls in that over 50% of such sources show 
signs of absorption in the range 0.7-0.9 keV, most likely 
due to oxygen edges, but relatively few possess conspicuous 




500 1000 1500 

H/S FWHM (km/sec) 



2000 



Figure 5. The relative luminosity of the soft excess compared to 
the hard power law (as given in Table 3) versus the H/3 line width. 
Objects with no evidence for absorption are marked by circles and 
those with normal warm absorbers by crosses. The sources ex- 
hibiting anomalous 1.1-1.4 keV absorption are marked by squares 
and Mkn 507 (with possible cold absorption) is marked with a 
filled circle. The arrows indicate upper limits on the strength of 
the soft excess. 



soft excesses, at least in the ASCA bandpass (e.g., Reynolds 
1997; George et al. 1998). 

A further interesting result illustrated by Fig. 6 
is the fact that all but one of the NLSls with 
H/3 FWHM < 1000 km/s show evidence for some sort 
of absorption, either by ionized material or (as in one 
case) by neutral gas. Conversely none of the NLSls with 
H/3 FWHM > 1000 km/s exhibit signs of hne-of-sight ab- 
sorption. One possibility is that this represents a very 
sharply tuned geometry dependence. For example, if the ab- 
sorbing material originates in an outflow from the disc, it 
would only be seen in low inclination systems (for a dis- 
cussion of the "pole-on" model for NLSls, see Boiler et al. 
1996 and references therein). It is unclear, however, why 
such an incredibly tight upper-limit should apply to the H/3 
line width in such circumstances. Clearly further study of 
this effect using both a larger sample of NLSls and a more 
sophisticated modelling of the putative absorption features 
is merited. 



5 CONCLUSIONS 

Our analysis of the ASCA X-ray spectral data for a sample of 
22 NLSls has revealed a number of important characteristics 
of this subclass of X-ray source. Absorption from neutral gas 
intrinsic to the host galaxy appears to be a very rare occur- 
rence in NLSl. On the contrary most NLSls exhibit a soft 
excess superposed on the hard power-law continuum. Mod- 
elled as a blackbody component, this soft excess often con- 
tains a significant fraction of the X-ray luminosity in these 
sources, even in the restricted ASCA bandpass. In the most 
extreme cases this presents an argument against an origin 
for the soft excess in the reprocessing of the hard contin- 
uum. The underlying power-law spectra span a wide range 
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of slopes (r = 1.6 ~ 2.5) with a mean spectral index of 2.1, 
which is slightly steeper than the norm for BLSls. Iron Ka 
emission lines are detected at a modestly significance level 
in roughly half of the objects observed, but unfortunately 
the line properties such as centroid energy, equivalent width 
and intrinsic width are in general rather poorly constrained 
by the ASCA data. 

The X-ray spectra of NLSls below 2 keV often exhibit 
additional features, which in the current paper we have in- 
terpreted as largely due to absorption by ionized material 
in the line of sight. Three of the NLSls show evidence for a 
absorption in the 0.7-0.9 keV range, probably arising from 
Ovil and Ovill edges, "warm absorber" features which are 
commonly observed in BLSls. However, six NLSls also show 
spectral features which we take to be absorption in the 
1.1-1.4 keV bandpass arising from resonance absorption in 
highly ionized material. A potentially very important result 
is that all objects showing signs of absorption by highly ion- 
ized gas lie in the lower half of the H/3 line width distribution 
{i.e. the absorbed sources have H/3 FWHM < 1000 km s"^). 
The combination of high effective area, good spectral reso- 
lution and extended bandwidth afforded by future missions 
such as XMM, AXAF and Astro-E, should allow unambigu- 
ous identification of the spectral features which appear to 
characterise NLSl galaxies and also confirm whether ion- 
ized absorption systems are preferentially observed in NLSls 
with extremely narrow H/3 line widths. 
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6 APPENDIX 

Some comments on the spectral fitting of individual 
objects 



Ton S180 

It is difficult to differentiate between absorption and emis- 
sion at ~1.2 keV in Ton S180. A single absorption edge 
at a rest energy of E= 1.09 ± 0.03 keV with an optical 
depth T = 0.17 ± 0.03, yields x^/^ = 1072/986. A sec- 
ond edge in the spectrum (at E= 1.42 ± 0.06 keV with 
r = 0.13±0.03) improves the fit further {x'/u = 1059/985). 
There is some evidence for an additional edge at ~0.73 keV 
iX^ /'^ = 1048/983) presumably corresponding to Ovil. How- 
ever SIS-1 seems to underestimate the flux compared to SIS- 
below 0.8 keV, so any features in the spectrum below 0.8 
keV need to considered with caution. 

A significant improvement over the soft excess fit is 
given when a broad Gaussian emission component is added 
to the model near 1 keV. The best-fit parameters are line 
EW= 20 ± 8 eV, energy E= 0.94 ± 0.03 keV and intrinsic 
width a = 0.062 ± 0.025 keV, for which xV^^ = 1055/985. 
A MEKAL model also gives an acceptable fit to the data, 
yielding a plasma temperature of kT= 0.89 ±0.11 keV 
ix^/u = 1050/986). 



IH 0707-495 

An absorption feature definitely gives the best fit in 
IH 0707-495. A broad absorption line improves the fit signif- 
icantly ix^/i' = 404/335), as per Table 4. An equally good 
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can be obtained, but for 1 additional free parameter, us- 
ing two edges instead of a single Gaussian; the first edge is 
at E=1.09±0.03 keV with t = 0.63±0.27 and the second at 
E=0.90±0.03 keV with r = 0.5 ±0.2. Leighly et al. (1997b) 
also find evidence for one or more absorption features near 
1 keV. The addition of a Gaussian emission line in the range 
0.6-3 keV, instead of an absorption feature, does not provide 
any improvement in the fit. 

NGC 4051 

Previous investigations of the X-ray spectrum of NGC 4051 
have suggested the existence of a warm absorber. There is 
no evidence for additional spectral features in the first ob- 
servation (taken during the PV phase of the mission) but 
the second (A02) observation, which has much better sig- 
nal/noise ratio, does show significant features in the residu- 
als between 0.7 and 0.9 keV. In this case we have modelled 
this apparent warm absorber with 2 absorption edges, for 
the first edge E= 0.73 ± 0.02 keV and r = 0.33 ± 0.05, for 
the second edge E=0.93 ± 0.03 keV and r = 0.19 ± 0.04. 
The resulting fit statistic is x^/i^ = 2010/1842, an improve- 
ment of Ax^ ~ 77 over the model without absorption. The 
results suggest a normal warm absorber in NGC 4051, with 
the edges probably originating from Ovil and Ovill. Our 
results are broadly consistent with those of Guainazzi et al. 
(1996), who consider the fitting of this observation in much 
greater detail. 

PG 1244 +026 

Emission and absorption components give comparable fits 
in PG 1244-1-026. The addition of a broad Gaussian absorp- 
tion line significantly improves the fit {x^ jv = 601/532), 
as per Table 4. A single edge at an energy of E=1.18 ± 
0.03 keV, gives a better fit {x jv = 591/533). Adding 
another edge further improves the fit; with two edges (at 
E=1.16±0.04 keV and ^=0.63^011 keV) the fit obtained is 
X^/v = 583/531. 

A broad Gaussian emission line at an energy of 
E=0.97±0.04 keV also significantly improves the fit jv = 
589/532), with an equivalent width of EW=36 ±11 eV. 
A slightly better fit (y^ /v = 586/533) is obtained with 
the MEKAL model, with a plasma temperature of kT=1.0 ± 
0.1 keV. Fiore et al. (1998) also find evidence for either an 
absorption feature at ~1.2 keV or an emission feature at 
0.9 keV. 



A broad Gaussian emission line at 0.74 ±0.05 keV gives 
a slightly worse fit {x^ /v = 470/423), as does the mekal 
model {x^ /v — 480/424) with a plasma temperature of 
kT=0.63 ± 0.05 keV. These results are broadly consistent 
with those of Leighly et al. (1997b). 



PG 1404+226 

Significant evidence is found for an absorption feature in 
PG 1404±226. The addition of a single edge to the best- 
fitting continuum model (PL±BB) yields x^/i^ = 172/165, 
an improvement of Ax^ = 29, with an edge energy of 
E= 1.07 ± 0.03 keV and optical depth r = 0.84 ± 0.28, con- 
sistent with the results of Leighly et al. (1997b). A slight im- 
provement (Ax^ =8.1) is also achieved by adding a second 
edge, although the edge energy of E= 2.4 ± 0.3 keV seems 
rather high. The 1 keV feature can also be modelled with 
a broad Gaussian absorption line, giving x^/i^ = 170/164 
with fit parameters as in Table 4. 

An emission model (Gaussian line) is unacceptable as 
the broad emission feature essentially models the continuum 
soft excess and gives an unacceptably large equivalent width. 



IRAS 17020+454 

The residuals from a simple power law plus Galactic absorp- 
tion fit show clear signs of absorption in IRAS 17020±454 
(see Figure 2). The addition of an edge at an energy of 
E=0.71 ± 0.02 keVwith a depth r = 0.91 ± 0.22, improves 
the fit significantly(x^/i^ = 979/971). A second edge at 
E=1.15 ± 0.04 keV with a depth r = 0.20 ± 0.08 further 
improves the fit jv ~ 965/969). A broad Gaussian ab- 
sorption line gives a comparable fit as detailed in Table 4. 

A narrow Gaussian emission line at an energy of 
E=1.06 ± 0.02 with an equivalent width of EW=23+7 eV 
gives a worse fit than an absorption line {x? /v = 1006/971). 
However, the inclusion of an edge at 0.72 ±0.02 keV, as well 
as the emission line, does lead to a net improvement in the 
fit {x^/u = 962/969). 

In agreement with Leighly et al. (1997a) we identify 
the 0.7 keV feature with an absorption edge from Ovil. 
Komossa & Bade (1998) present a ROSAT spectrum of 
IRAS 17020±454 and also find evidence for a warm ab- 
sorber, albeit a dusty one. 



IRAS 13224-3809 

Absorption components give a slightly better fit than emis- 
sion components. The addition of a single absorption edge to 
the model improves the fit significantly {x^/v = 478/424). 
The best fit edge energy is E=l.ll ± 0.04 keV. A second 
edge improves the fit further (x^ jv = 467/422), with edge 
energies at E=1.00 ± 0.04 and 1.18 ± 0.03 keV and optical 
depths of r = 0.73 and 1.15 respectively. This fit also al- 
ters the continuum parameters slightly: F = 1.84 ±0.07 and 
kT=163 ± 9 eV. A broad Gaussian absorption line gives a 
slightly better fit {x^ /y ~ 459/423), with parameters as in 
Table 4. 



IRAS 2OI8I-224 

There may be an emission or absorption feature in the spec- 
trum of IRAS 20181-224, but the improvement in the fit 
upon adding either is rather small. A broad absorption line 
gives a slight improvement in the fit, as noted in Table 4, 
but adding one or two edges in the range 0.6-3 keV gives a 
significantly a worse fit. 

A narrow Gaussian emission line at 1.11 ± 0.05 keV 
also slightly improves the fit {y^ /u = 438/387). The mekal 
model gives a similar result with a plasma temperature of 
kT=1.7to:3 keV. 
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Ark 564 

Comparable fits are obtained witli eitlier a broad emission 
feature at ~1 keV or two absorption components. A sin- 
gle absorption edge at E=1.23 ± 0.02 keV gives x^/'^ ~ 
1862/1565 whereas two edges improves the fit to X^/*^ ~ 
1819/1563 (the edge energies being E=1.14 ± 0.03 keV and 
1.37±0.04 keV with r = 0.18±0.04 and 0.15±0.03, respec- 
tively). This is consistent with the Brandt et al. (1994) mod- 
elling of the ROSAT spectrum of Ark 564 which employed 
a deep absorption edge at 1.15 keV. The broad Gaussian 
absorption line model detailed in Table 4 is clearly not the 
best representation of the data {"x^ /v = 1854/1564). 

A broad emission line provides a slight improvement 
over two edges {x^/v = 1815/1564) with E=0.99±0.01 keV, 
(J = 0.10 ± 0.03 keV and EW=39l| eV. The mekal model 
gives a worse fit [x^ /v = 1830/1565), with a plasma tem- 
perature of kT=1.10l^;o^ keV. 

The SIS and GIS detectors systematically diverge be- 
low 0.9 keV. There are also signs that the two SIS detectors 
are diverging at ^0.6 keV with SIS-1 seeming to underesti- 
mate the flux compared with SIS-0. These problems make 
it difficult to establish the nature of the soft X-ray spectral 
features in Ark 564. 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science ETIeX style file. 
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